. KPP Front Speeds in Random Shears and 

the Parabolic Anderson Problem 

; Jack Xin 

Q ; Department of Mathematics k TIC AM 



On 



X 

5-S 



University of Texas at Austin 
Austin, TX 78712, USA. 



> 

in 
o 

<N 

(N 
O 

Q^c Abstract 

We study the asymptotics of front speeds of the reaction-diffusion 
equations with Kolmogorov-Petrovsky-Piskunov (KPP) nonlinearity 
and zero mean stationary ergodic Gaussian shear advection on the 
entire plane. By exploiting connections of KPP front speeds with the 
almost sure Lyapunov exponents of the parabolic Anderson problem, 
and with the homogenized Hamiltonians of Hamilton-Jacobi equa- 
tions, we show that front speeds enhancement is quadratic in the small 
root mean square (rms) amplitudes of white in time zero mean Gaus- 
sian shears, and it grows at the order of the large rms amplitudes. 
However, front speeds diverge logarithmically if the shears are time 
independent zero mean stationary ergodic Gaussian fields. 



1 Introduction 



Front propagation in random media appears in many scientific areas, and 
analysis of prototype models has been an efficient way to improve our under- 



standing of fronts with complex structures ||21|| . We shall consider reaction- 
diffusion (R-D) fronts moving inside a two dimensional random shear field 
modeled by: 

ut = - AyU + at- Vu + f(u), y = ( yi ,y 2 ) G R\ (1.1) 

where £ = (0, t))), £ a mean zero stationary ergodic Gaussian field, 
a > a parameter measuring the shear strength (root mean square amplitude 
for short), A y = dy +dy 2y2 , f(u) = u(l— u), known as the KPP nonlinearity. 
We are interested in fronts along the y 2 direction, and how the multi-scales 
in the shear field £ influence the speed of propagation, a quantity at large 
scale (large time). 

If the shear £ is deterministic, e.g. periodic in space (and time), the 
asymptotic behavior of front speed is known for many forms of reaction 
nonlinearity f(u). For small a and zero mean £: 

c* = c (l + aa 2 + h.o.t), (1.2) 

where a is a positive constant depending only on £, Co is the pure R-D front 
speed. The constant a is universal among bistable nonlinearity (f(u) = 
u(l — u){u — /i), fi G (0,1/2)), ignition combustion nonlinearity (f(u) = 
0,u G (0,0), 9 G (0,1), ignition level, f(u) > 0, u G (0,1)), and others 
(f(u) = u m (l — u), m > 2), except KPP (f(u) = u(l — u)) where a takes 
another positive value. See calculations using perturbation and variational 



methods [ 16| , [ij], in the spatially periodic case. The robustness of ( fL2 ) in 



the spatially-temporally periodic case is studied recently in [H]. Except for 



the KPP nonlinearity, where linearization is possible to obtain front speeds, 
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one needs to work with an exact traveling wave solution and the resulting 
equation it satisfies in order to find the constant a. 
For large a, on the other hand: 

c* = 0(a), (1.3) 

if/(w) > 0, see !, 0, 0, 0, |g, and existence of limit c* ja holds for the 
KPP [|[. The asymptotics Q1.3Q also hold in the spatially-temporally periodic 
shears [ 14] , and is numerically observed even for bistable nonlinearity |jT4] . 



A natural question is whether these results remain valid for fronts in ran- 
dom shears. We approach this issue here for the KPP nonlinearity and in 
the entire space R 2 . We show that the KPP front speed formula is connected 
with the almost sure Lyapunov exponent of the so called parabolic Ander- 
son problem, namely the exponential growth rate of the semigroup e tL with 
L- — A ■ +V(x,t)-, V(x,t) a Gaussian process with mean zero. There is a 
rich literature on almost sure Lyapunov exponent [f§], ||, ||. The logarithm 
of the linear parabolic equation satisfies the stochastic viscous KPZ equation 
in surface growth (or a quadratic Hamilton- Jacobi equation with random 
potential and viscosity), and the growth rate is also the effective Hamilto- 
nian in the large space large time (homogenization) limit. Combining these 
two connections, we are able to simplify the front speed calculations. The 
main finding is: if £ = £{yi,t), white in time, mean zero stationary ergodic 
Gaussian process, then almost surely front speed in random shears is a de- 
terministic constant c* with asymptotics: 

c* = c (l + ^T{0)a 2 + h.o.t), <7«1, (1.4) 

T(-) the spatial covariance function of £, T(0) > 0; and: 

c* = 0(a), ct>1., (1.5) 

However, if £ = £(?/i), time independent mean zero stationary ergodic 
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Gaussian process, then almost surely, the front velocity diverges like: 



c ~ 



c + aJ2T(0) log t, 



(1.6) 



due to the almost sure growth of running maxima of Gaussian process. The 
divergence can be faster for other types of noisy processes, for example the 
Poisson shot noise. 

(1.5) and ( |1.6|) , we see that the fast switching of direc- 



Comparing (|T 



tion in time helps to slow down fronts, and cures the divergence in (|1.6|) . The 
slow down of front speeds due to fast direction changing in time also occurs 



in the deterministic setting for all the above nonlinearities JL4| . As the front 
speed asymptotics in the deterministic cases (|1.2|) - (|L3|) are not sensitive to 
nonlinearities, we expect that the analysis on random KPP front speeds, 
(|1.4| )- (|1.6|) , will have implications for other types of f(u). Numerical work 



along this line is in progress and will be reported elsewhere ITS 



The rest of the paper is organized as follows. In section 2, we derive the 
front speed asymptotic formula ( |1.4j) , (|1.6|), and bound (|1.5|), using properties 
of homogenized Hamiltonian and asymptotics of Lyapunov exponents. In 
section 3, we draw conclusions and discuss future works. 



2 KPP Front Speeds in Random Shear Flows 

The KPP (minimal) front speeds in random media have a variational char- 
acterization, ||, jnj, [21|, and references therein, only /'(0) enters from 
nonlinearity in an otherwise linear problem. 

Let £ = £(yi,t,u), a mean zero stationary spatial-temporal Gaussian 
field, with covariance function: 

E(£(yu t) ^y'v *')) = r(|t - t'\, \y x -y[\), (2.1) 
V y, y' G R 1 , t, t' G R 1 , where T is a positive definite continuous function 
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with decay at infinity. The KPP front speed formula is: 

c* = inf H ^ X *\ (2.2) 
A 2 >o A 2 

where: 

H(X) = /'(0) + lim r 1 log E[exp{\ ■ Y t }\, (2.3) 

t — >oo 

and Y t is a diffusion process obeying the Ito equation: 

dY s = dW s + (0, a£(y M , t - s)) ds, Y = 0, (2.4) 

W s is the standard two dimensional Wiener process. 
It follows from ( ^4|) that: 

dF 2iS = tW 2 , a + a£(W M ,t-s)ds. 

So: 

The second factor = exp{^|-^}. The first factor = u(0, t), with it = u(x,t) 
the solution of: 

u t = -u xx + \ 2 cr£(x,t)u, xeR 1 , 

u(x,0) = e XlX . (2.6) 

The initial value problem is invariant in the sense of distribution if (Ai, A2) — > 
— (Ai, A2), x — > —a;, implying that the limit: lim^oo logw(0, £) = 7(A) is 
an even function in A, if it exists and is nonrandom. 

The function v = \ogu(x,t) satisfies the viscous Hamilton- Jacobi equa- 
tion: 

v t = -v xx + -vl + \ 2 a£(x,t), (2.7) 
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with linear initial data v(x,0) = X±x. The limit lim^oo t 1 v agrees with 
the homogenized Hamiltonian of the related homogenization problem, see 



21|] for more exposition of this connection. A useful consequence is that 
7 is a convex function in Ai being the homogenization limit of a quadratic 
Hamiltonian with oscillating potential (for fixed A 2 ), see |17| , fll8l , jl9| for 
analysis of convex Hamiltonian and bounded random potentials. 

Evenness and convexity of 7 in Ai implies that 7(0, A2) < 7(Ai, A2), and: 

A 2 >o A2 

As Ai = 0, the problem reduces to the parabolic Anderson model, studied by 
Carmona, Molchanov, and coworkers, |4|], 0, ||. The limit, 

lim t~ l logw(0,t) =7, 

if it exists and is nonrandom, is called the almost sure Lyapunov exponent. 

A case where the 7 limit exists as a finite nonrandom number is when 
£ is white in time, or T{t,yi) = 5{t)T (\yi\). The spatially discrete case on 
lattice Z d was completely analyzed in |4[], ||, where the Laplacian operator 
A x v is replaced by J2\ x '~x\=i v(x') — v (x). The main result is: Let v be the 
solution of 

v t = kAv + £(t, x) v, x G Z d , k>0, (2.9) 

with product in the Stratonovich sense, then the almost sure Lyapunov ex- 
ponent j(k) exists as nonrandom number and obeys (60 > 0, d — 1, 2): 

p— ' K<<1 ' ( 2 - 10 ) 

— log K 

7(«) ~ ^ «>1- (2.11) 

In three and higher dimensions, 7(«) = T(0)/2 if k exceeds a critical level k ct . 
This shows that the larger the diffusion effect, the more 7 tends to be constant 
T(0)/2. Moreover, j(k) is a monotone increasing continuous function in k [ffl. 
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The continuum case is more involved, Cranston and Montford just proved M 
that j(k) ~ cik p , p e (0, 1/2), ci > 0, k < 1, and that 7(«) > for all k > 0. 
The large k regime ( |2.11| ) is expected to be the same, also j(k) is monotone 
in k. If the continuum A is approximated by standard central differencing, 
we see that the continuum A is close to a discrete one with larger diffusion. 
The results can be adapted to equation (|2.6| ) with Ai = 0. First ^£(i/a, y) = 
£(t, y) in law, for any a > 0. Let t — t' / a, a — (<tA 2 ) 2 , equation fl2.6|) becomes: 



M ( , = ^A« + ^',i)tt, (2.12) 
It follows that lim^oo u(t, 0)/t — 7*(crA 2 ) such that 

7>A 2 ) ~ r( ° )( 2 aA2)2 ; aA 2 «l, (2.13) 

7 >A 2 ) ~ b h ^^5) ( ffA 2) 2 ' ffA 2 » 1, (2.14) 

where h = h(x) = x p , p G (0, 1/2). Now we minimize H(0, A 2 )/A 2 . 
For a -C 1, and A 2 ~ 0(1), where the minimizer is expected: 

H(0, A 2 )/A 2 ~ /'(0)/A 2 + A 2 /2 + T ^f X2 + h . . tj 



giving minimal point: 



and minimum: 



A, 



2f(0) 

\ 1 + r(o)a 2 
1. 



2f (0) (1 + |r(0)a 2 ) + h.o.t. (2.15) 

This is the stochastic analogue of quadratic speed enhancement. 
For a 1, and A 2 <r ^> 1, we have: 

H(0, A 2 )/A 2 ~ b h J ^ 2 A 2 + f(0)/A 2 + A 2 /2, 
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giving an upper bound with A2 = — , ctq large enough constant so that 
a > d > 1: 

c* = inf #(0,A 2 )/A 2 < 0(al- 2p a), (2.16) 
A 2 >o 

which shows that the upper bound on the speed is linear in a ^> 1. 

To establish a similar lower bound, let us note that j(k) is bounded away 
from zero by an amount 70 if k > k q , and asymptotics j(k) ~ C\K V holds, for 
some p G (0, 1/2), if k < kq. 

Now we estimate c* from below. If 2(<rA 2 ) 2 < Kq 1 , then: 

H(0, A 2 )/A 2 > 7o <r 2 A 2 + f(0)/A 2 + A 2 /2 = F(A 2 ). (2.17) 
The absolute minimal point A 2jm of F satisfies: 

l/2-/'(0)/A 2 + 7o^ 2 = 0, 

or: 

x ^=0S~ oia ^ a>>i - (2 - i8) 

So A 2 , m( 7 = 0(1), and H(0, A 2 )/A 2 > F(A 2 ) > 0(a). 
On the other hand, if 2(<rA 2 ) 2 > Kq 1 , 

H(0, A 2 )/A 2 > Cl 2- p a 2 - 2p \\- 2p + f(0)/A 2 + A 2 /2 = G(X 2 ). (2.19) 

The absolute minimal point \ 2 ,m of G obeys: 

ci2-*(l - 2p)a 2 ~ 2p \f p + 1/2 = /'(0)/A 2 , 

implying: 

Cl 2"'(l - 2p)a 2 - 2p \~% < /'(0)/A| M , 

or: 

A 2 ,M<C(p, Cl ,f(0))/a, (2.20) 
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for some positive constant C(p, ci, f'(0))- It follows from ( |2.20| ) and ( f2.19| ) 
that: 

#(0,A 2 )/A 2 > f(0)/A 2 , M > C"(p, Cl ,/'(0))<7. 

The above two lower bounds combine to give c* > 0(a), and together with 
upper bound ( [2.161 ), we have the front speed bound c* = 0(a) for large a. 

The 7 diverges logarithmically when £ = is a spatial Gaussian process 
|J. The leading order exponential growth rate of u depends on the fact that: 



sup i(x) ~ V 2r (°) lo S*> 1 °°' a - s -' ( 2 - 21 ) 

|sj|<* 

and so: 

lim = log u(t,x) = a|A 2 K/2r(0), a.s. (2.22) 

implying front speed divergence in time. In other words, at large t: 

A 2 A 2 z 

or: 

c * = in f _ Cq + av / 2 r(0)logt. (2.23) 

A 2 >o A 2 

In view of ( |2.21| ), such divergence can be avoided if the front is restricted 
inside an infinite channel of finite cross section. However, the almost sure 
constant speed ceases to exist, as on a finite interval £ process loses ergodicity, 
different realizations can be either more positive or negative, so ensemble 
averaged speed becomes an alternative measure of front speed. A numerical 
study of random fronts in channels will be reported elsewhere [ff5 . 



3 Conclusions 

We found that KPP front speeds are finite and enhanced through white in 
time stationary Gaussian random shears. The enhancement is quadratic in 
the root mean square (rms) of small shear amplitudes, and grows at the 
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order of the large rms amplitudes. The front speeds diverge logarithmically 
in time independent mean zero Gaussian shears. The findings will help to 
investigate similar front behaviors under other forms of nonlinearities, and 
serve as a useful guide for numerical works on fronts through random shears 
in channels of large but finite cross sections. 
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